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We present the results of searches for decays of B mesons to final states with a bi meson and 
a charged pion or kaon. The data, collected with the BABAR detector at the Stanford Linear 
Accelerator Center, represent 382 million BB pairs produced in e + e~~ annihilation. The results for 
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the branching fractions are, in units of 10" 6 , B(B+ -> 6?vr + ) = 6.7±1.7±1.0 (A.Oct), B(B + -> b°K + ) 
= 9.1±1.7±1.0 (5.3cr), B(B° -» fcf tt ± ) = 10.9±1.2±0.9 (8.9<r), and B(B° -> 67/^+) = 7.4±1.0±1.0 
(6.1cr), with the assumption that B(bi — > lj7t) = 1. We also measure charge and flavor asymmetries 
Ach(B+ -> 6?7r+) = 0.05±0.16±0.02, A ch (B+ -> b^K+) = -0.46±0.20±0.02, A c h(B° -> bfir ± ) = 
-0.05 ±0.10 ±0.02, C(B° -> 6f7r ± ) = -0.22 ±0.23 ±0.05, AC(B° -» 6f tt*) = -1.04 ±0.23 ±0.08, 
and Ach(B° -» 67if + ) = -0.07 ± 0.12 ± 0.02. The first error quoted is statistical, the second 
systematic, and for the branching fractions, the significance is given in parentheses. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



Recent searches for decays of B mesons to final 
states with an axial-vector meson and a pion have re- 
vealed modes with rather large branching fractions, e.g., 
B(B° -► afir*) = (33.2 ± 3.8 ± 3.0) x 10" 6 [j]. Here we 
search for related modes with a 6? or a meson plus 
a 7T+ or K+ 0, in a sample of (381.8 ± 4.2) x 10 6 BB 
pairs produced by e + e~ annihilation at the T(AS) res- 
onance (center-of-mass energy y/s — 10.58 GeV). The 
integrated luminosity is 346 fb _1 . 

The mass and width of the b t are 1229. 5±3. 2 MeV and 
142 ± 9 MeV, respectively, and the dominant decay is to 
ujir Q. In the quark model the b\ is the I G — 1 + mem- 
ber of the J PC = 1 H , X P\ nonet, whereas the ai is the 
I G = 1" state in the J PC = 1++, 3 Pi nonet. The avail- 
able theoretical estimates of the branching fractions of 
B mesons to &i7r and b\K come from calculations based 
on naive factorization 0)0), and on QCD factorization 
Q. The latter incorporate light-cone distribution ampli- 
tudes evaluated from QCD sum rules. Expected branch- 
ing fractions lie in the range 5-10 x 10 -6 [6j; estimates as 
large as 26xl0~ 6 are found in the calculations of and 
40xl0~ 6 in those of @|. 

The four modes B+ -> 6?7r+, B+ -> b\K + , B° -> 
b^TT + , and B° — > b^K + can be mediated by external 
tree amplitudes in which the weak current produces the 
pion (kaon) with a Cabibbo-favored (suppressed) cou- 
pling. Alternatively, a "penguin" loop amplitude is fa- 
vored for the kaon modes, and suppressed for the pion 
modes. The fifth mode, B° — > b^n" , requires a cou- 
pling of the current to the bf , which is forbidden for this 
G = +1 state 0, leading to the expectation B(B° — > 

6+7T-) < B(B° -> b^TT+). 

Direct CP violation would be indicated by a non-zero 
value of the asymmetry A ch = (r - - r+)/(r - in 
the rates T ± (B ± — > F^) for decay of a charged B meson, 
or T + (B° — > b^K + ) and its charge conjugate. For the 
decays B° — > bfir^ 1 we define A c h and two additional 
asymmetries C and AC through 

r,,/ = J(r + r) (1 ± [i + / (c + q AC)\ , (1) 

where the signal B meson flavor / = ±1 for B°, — 1 for 
B°, and q is the sign of charge of the b\. To measure C 
and AC we use the flavor r\ (±1 for B° and —1 for B°) 
of the second meson Stag produced in T(4S I ) decay Q. 



The yields are given by 

Y m = -Y s (l + qA ch ){l-riAw + r)n(l-2w) (2) 

- 77(1 - 2 Xd ) [l-2w + fi(r) - Aw)} (C + qAC)} , 

where Y$ is the total signal yield, \d = 0.188 ± 0.003 
the time-integrated mixing probability Q , w the mistag 
fraction, and Aw and \x the B — B differences in the 
mistag rate and tagging efficiency, respectively. 

The data were collected with the BABAR detector 9] at 
the PEP-II asymmetric e + e~ collider located at the Stan- 
ford Linear Accelerator Center. Charged particles from 
the e + e~ interactions are detected, and their momenta 
measured, by a combination of five layers of double-sided 
silicon microstrip detectors and a 40-layer drift chamber, 
both operating in the 1.5 T magnetic field of a supercon- 
ducting solenoid. Photons and electrons are identified 
with a CsI(Tl) electromagnetic calorimeter (EMC). Fur- 
ther charged particle identification (PID) is provided by 
the average energy loss (dE/dx) in the tracking devices 
and by an internally reflecting ring imaging Cherenkov 
detector (DIRC) covering the central region. A detailed 
Monte Carlo program (MC) is used to simulate the B 
production and decay sequences, and the detector re- 
sponse [lOj. 

The bi candidates are reconstructed through the de- 
cay sequence b\ — > lot:, ui — > 7r + 7r~7r°, and tt° — > 77. 
The invariant mass of the photon pair is required to lie 
between 120 and 150 MeV, i.e., within about two stan- 
dard deviations of the nominal mass [3j . For the b\ and 
ui whose masses are observables in the maximum likeli- 
hood (ML) fit described below, we accept a range that 
includes wider sidebands (see Fig.Q]). Secondary charged 
pions in b± and ui candidates are rejected if classified as 
protons, kaons, or electrons by their DIRC, dE/dx, and 
EMC PID signatures. For the primary pion (kaon) from 
the J3-meson decay we define the PID variable S n (Sk) 
as the number of standard deviations between the mea- 
sured DIRC Cherenkov angle and that expected for a 
pion (kaon), requiring — 2 < < 5 (— 5< Sk < 2). 

We reconstruct the S-meson candidate by combining 
the 4- momenta of a pair of daughter mesons, using a 
fit that constrains all particles to a common vertex and 
the 7T° mass to its nominal value. From the kinemat- 
ics of T(AS) decay we determine the energy-substituted 

mass mEs = \ \ s ~ Ps and energy difference AE = 
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Eb — j\/s, where (Eb,Pb) is the B-meson 4-momentum 
vector, and all values are expressed in the T(45) rest 
frame. The resolution in toes is 2.4 — 2.7 MeV and in 
AE is 25-32 MeV, depending on the decay mode. We 
require 5.25 GeV < to es < 5.29 GeV and -0.13 GeV < 
AE < AE max , with AE max = 0.1 (0.13) GeV for b\ (&+), 
where the tighter restriction serves to limit the number 
of combinatorial candidates per event. 

We also impose restrictions on resonance decay an- 
gles to exclude the most asymmetric decays where soft- 
particle backgrounds accumulate and the acceptance 
changes rapidly. We require cos 6^ < 1.1 — 0.5 1 cos 0^1, 
where 9i >1 is the angle between the momenta of the pion 
from b\ — > ton and its parent B meson, measured in the b\ 
rest frame, and 9^ is the angle between the normal to the 
lu — ► 3-7T decay plane and the momentum of its parent b\ , 
measured in the u) rest frame. Backgrounds arise primar- 
ily from random combinations of particles in continuum 
e + e~ — > qq events (q = u, d, s, c). We reduce these with 
a requirement on the angle 9t between the thrust axis 
of the B candidate in the T(4S) frame and that of the 
rest of the charged tracks and neutral calorimeter clus- 
ters in the event. The distribution is sharply peaked near 
| cos#t| = 1 for qq jet pairs, and nearly uniform for B- 
meson decays. The requirement, which optimizes the ex- 
pected signal yield relative to its background-dominated 
statistical error, is |cos#t| < 0.7. The average number 
of candidates found per selected event is in the range 1.3 
to 1.4 (1.4 to 1.6 in signal MC), depending on the final 
state. We choose the candidate with urn invariant mass 
closest to the nominal value of the b\ mass 0] . In the ML 
fit we discriminate further against qq background with 
a Fisher discriminant T that combines several variables 
which characterize the energy flow in the event [llj . It 
provides about one standard deviation of separation be- 
tween B decay events and qq background. 

We obtain yields for each channel from an extended 
ML fit with the input obscrvables AE, imes, T , and the 
resonance masses rrib 1 and to w . The selected data sample 
sizes are given in Table HI Besides the signal events these 
samples contain qq (dominant) and BB with b — > c 
combinatorial background, and a fraction of cross feed 
from other charmless BB modes, which we estimate from 
the simulation to be (0.5-0.8)%. The last include non- 
resonant u!tt(tt,K), and modes that have final states dif- 
ferent from the signal, but with similar kinematics so 
that broad peaks near those of the signal appear in some 
observables, requiring a separate component in the prob- 
ability density function (PDF). The likelihood function 
is 

N 

C = cxp (- ]T Y jm ) H £ Y im x (3) 

Vj (m ES 1 )Vj (P )Vi (ae* )Vi K x )v 3 K ) , 

where N is the number of events in the sample, and for 



each component j (signal, combinatorial background, or 
charmless BB cross feed), Yj iqn is the yield of events 
(Eq. [2|) and Vj(x l ) the PDF for observable x in event i. 
The signal component is further separated into two com- 
ponents (with proportions fixed in the fit for each mode) 
representing the correctly and incorrectly reconstructed 
candidates in events with true signal, as determined with 
MC. The factored form of the PDF indicated in Eq.[3]is 
a good approximation, particularly for the combinatorial 
qq component, since we find correlations among observ- 
ables in the data (which are mostly qq background) to be 
small. The effects of this approximation are determined 
in simulation and included in the bias corrections and 
systematic errors discussed below. 

We determine the PDFs for the signal and BB back- 
ground components from fits to MC samples. We cali- 
brate the resolutions in AE and toes with large data con- 
trol samples of B decays to charmed final states of simi- 
lar topology (e.g. B — > D(Kirir)ir). We develop PDFs for 
the combinatorial background with fits to the data from 
which the signal region (5.27 GeV < toes < 5.29 GeV 
and \AE\ < 0.1 GeV) has been excluded. 

The functions Vj are constructed as linear combi- 
nations of Gaussian and polynomial functions, or in 
the case of toes for qq background the threshold func- 
tion xyl — x 2 exp [— £(1 — a; 2 )] , with argument x = 
2toes/v^ an d parameter £. These functions are dis- 
cussed in more detail in [llj], and are illustrated in Fig.Q] 

We allow the parameters most important for the deter- 
mination of the combinatorial background PDFs to vary 
in the fit, along with the yields for all components, and 
the signal and qq background asymmetries. Specifically, 
the free background parameters are: £ for toes, linear 
and quadratic coefficients for AE, and the mean, width, 
and width difference and polynomial fraction parameters 
for T. 

We validate the fitting procedure by applying it to en- 
sembles of simulated experiments with the qq component 
drawn from the PDF, into which we have embedded the 
expected number of signal and BB background events 
randomly extracted from the fully simulated MC sam- 
ples. Biases obtained by this procedure with inputs that 
reproduce the yields found in the data are reported, along 
with the signal yields, in Table [U 

In Fig. Q] we show the projections of the PDF and data 
for each fit. The data plotted are subsamples enriched 
in signal with a threshold requirement on the ratio of 
signal to total likelihood (computed without the plotted 
variable) that retains (29-53)% of the signal, depending 
on the mode. 

We compute the branching fraction by subtracting the 
fit bias from the measured yield, and dividing the result 
by the efficiency times B{u -> tt+tt-tt ) = 89.1±0.7% 0, 
and by the number of produced BB pairs. We assume 
r(T(45) -> B+B-)/T(T(AS) -> B°B°) = 1, consistent 
with measurements 0]. The results are given in Table U 
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TABLE I: Number of events N in the sample, fitted signal yield Is, and measured bias (to be subtracted from Yg) in 
events (ev.), detection efficiency e, significance S (with systematic uncertainties included), and branching fraction and charge 
asymmetry with statistical and systematic error. 



Mode N (ev.) Y s (ev.) Bias (ev.) e (%) S (cr) B (10 -6 ) 



<Ach 



b?7T + 


32f76 


178 +3 3 9 7 


26±14 


6.78 


4.0 


6.7 ± 1.7 ± 1.0 


0.05 ±0.16 ±0.02 


blK+ 


f8036 




24±12 


6.73 


5.3 


9.1 ± 1.7 ± 1.0 


-0.46 ±0.20 ±0.02 


bfn ± 


3690f 


387I« 


34±17 


9.54 


8.9 


10.9 ±1.2 ±0.9 


-0.05 ±0.10 ±0.02 


b^K + 


17497 


267lg 


32±16 


9.43 


6.1 


7.4 ± 1.0 ± 1.0 


-0.07 ±0.12 ±0.02 



. ^ 20+; 
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FIG. 1: Distributions for signal-enhanced subsets of the data projected onto the fit observables for the decays: (a-e) B + — > fe?7r + , 
(f-j) B + — > b\K + , (k-o) B° — > bf-n^, and (p-t) B° — > b± K + . The solid line represents the result of the fit, and the dashed line 
the background contribution. 



along with the significance, computed as the square root 
of the difference between the value of — 21n£ (with ad- 
ditive systematic uncertainties included) for zero signal 
and the value at its minimum. 

Systematic uncertainties on the branching fractions 
arise from the PDFs, BB backgrounds, fit bias, and effi- 
ciency. PDF uncertainties not already accounted for by 
free parameters in the fit are estimated from the consis- 
tency of fits to MC and data in control modes. Varying 
the signal-PDF parameters within these errors, we esti- 
mate yield uncertainties of (2.4-3.3)%, depending on the 
mode. The uncertainty from fit bias (Table[l| includes its 
statistical uncertainty from the simulated experiments, 
and half of the correction itself, added in quadrature. For 
the BB backgrounds we vary the fixed fit component by 
100% and include in quadrature a term derived from MC 
studies of the inclusion of a b — ► c component with the 
dominant qq background. Uncertainties in our knowledge 
of the efficiency include 0.5% x N t and 1.5% x N~, where 



N t and iV 7 are the numbers of tracks and photons, re- 
spectively, in the B candidate. The uncertainties in the 
efficiency from the event selection are below 0.5%. 

We study asymmetries from the track reconstruction 
(found negligible), and from imperfect modeling of the 
interactions with material in the detector, by measur- 
ing the asymmetries in the qq background in the data 
and control samples mentioned previously, in compari- 
son with MC [l2| . We apply corrections, and assign sys- 
tematic errors, to A c h equal to —0.010 ±0.005 for modes 
with a primary kaon and 0.000 ± 0.005 for those with a 
primary pion. The leading systematic errors on C and 
AC come from the fit bias. 

With the assumption that B(bi — > lott) = 1, we obtain 
for the branching fractions: 



B(B+ 
B{B+ 
B{B° 
B(B°- 



> b\K+) 
• bf^) 
KK+) 



(6.7 ± 1.7 ± 1.0) x 10~ 6 
(9.1 ±1.7 ±1.0) x 10~ 6 
(10.9 ± 1.2 ±0.9) x 10~ 6 
(7.4 ± 1.0 ± 1.0) x 10~ 6 . 
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For the asymmetries we find 



A ch (B+ 


-+&°7T+) = 


0.05 ±0.16 ±0.02 


A ch (B+ - 


^b\K+) = 


-0.46 ±0.20 ±0.02 


A ch (B a - 




-0.05 ±0.10 ±0.02 


C(B° - 


= 


-0.22 ±0.23 ±0.05 


AC(B° - 


-+&?7r±) = 


-1.04 ±0.23 ±0.08 


A ch (B° - 


+ KK+) = 


-0.07 ±0.12 ±0.02. 



The first error quoted is statistical and the second sys- 
tematic. The QCD factorization estimates Q for the 
branching fractions and charge asymmetries agree with 
these measurements within experimental and theoreti- 
cal errors. The authors of [6j note that the observa- 
tion B(B+ -► b° 1 K+)/B(B° -> b^K+) > 0.5, if con- 
firmed with higher precision, would indicate the presence 
of a weak annihilation contribution to these modes. The 
value of the CP-conserving AC near —1 for B° — > bf^ 
agrees with the expected suppression of B° — > bfir~; our 
results imply the ratio T(B° -> &+7r-)/r( J B° -> 6] F 7r ± ) = 
— 0.01±0.12. We find no evidence for direct CP violation 
in these decays. 
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